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ABSTRACT 

Tills  investigation  is  for  the  study  of  the  electronic 
structure  of  actinide  complexes.   Accurate  Hartree-Fock  self- 
Gonslstent-f ield  computations  have  been  performed  in  the 
recently  developed  X  -Multiple-Scattering  approximation.   The 
goal  of  the  investigation  is  to  interpret  the  optical  and 
magnetic  properties  of  the  actinide  compounds  and  to  estimate 
the  degree  to  which  covalency  must  be  taken  into  account  in 
order  to  understand  the  role  of  the  5f  electrons  in  both 
molecular  and  metallic  bonding. 

The  compounds  to  be  investigated  will  include  the  actinyl 
ions  as  well  as  the  hexafluoride  and  oxide  complexes  of  Uranium, 
Neptunium,  and  Plutonium.   In  addition,  it  is  proposed  to  inves- 
tigate small  metallic  clusters  which  reflect  the  complex  geometry 
of  the  unit  cell  found  in  the  metallic  phase  of  the  actinide 
elements . 


I.   Introduction 

This  report  outlines  some  of  the  preliminary  progress  made 
In  the  study  of  the  electronic  structure  of  actlnlde  complexes 
In  the  Hartree-Fock  self-consistent-field  approximation.   Our 
goal  is  to  interpret  the  optical  and  magnetic  properties  of  these 
compounds  and  to  estimate  the  degree  to  which  covalency  must  be 
taken  into  account  in  order  to  understand  the  role  of  the  5f 
electrons  in  both  molecular  and  metallic  bonding.   The  specific 
compounds  to  be  investigated  will  include  both  the  hexafluoride 
and  oxide  complexes  as  well  as  small  metallic  clusters  which 
reflect  the  complex  geometry  of  the  unit  cell  found  in  the 
metallic  phase  of  these  elements. 

Historically  the  spectra  of  actlnlde  compounds  have  been 
discussed  in  terms  of  crystal  field  theory  in  which  the  optical 
spectra  is  interpreted  in  terms  of  the  Slater-Condon  parameters 
(arising  from  electronic  repulsion),  the  spin  orbit  coupling 
parameter,  and  the  crystal  field  splitting  parameters  (due  to 
the  electrostatic  field  of  the  llgands ) .   On  comparison  with 
experiment  one  finds  both  the  Slater-Condon  parameters  and  the 
spin  orbit  parameters  are  reduced  from  the  free  atom  values. 
This  change  from  the  free  atom  parameters  is  interpreted  to 
arise  from  the  covalent  mixing  of  ligand  orbitals  into  the 
metal  d  and  f  functions.   The  covalency  effects  in  actlnlde 
complexes  are  generally  thought  to  be  intermediate  between  that 

of  the  well  shielded  ^f   electrons  of  the  rare  earths  and  the 

1 ,2 
strong  covalent  bonding  found  in  many  transition  metal  ions. 


Indeed,  recently,  Hlll-^  has  argued  that  the  5f  electrons  are  often 
better  compared  to  the  3d  electrons  of  the  first  transition  series 

than  to  the  ^f  electrons  of  the  rare-earth  series.   The  results  of 

4 
recent  band  structure  calculations  by  Kmetko   and  Koelling  and 

5 
Freeman  would  seem  to  give  weight  to  this  hypothesis. 

In  an  attempt  to  understand  covalency  effects  in  transition 
metal  ions  we  have  recently  carried  out  Hartree-Fock  self- 
consistent-field  calculations  on  the  NiFg    cluster.    The 
molecular  orbitals  of  the  complex  were  expanded  in  a  basis  of 
Gaussian  functions  which  included  sufficient  flexibility  to 
allow  for  orbital  relaxation  in  the  field  of  the  ligands.   All  86 
electrons  were  included  and  all  multicenter  integrals  were 
accurately  calculated.   Even  for  this  ionic  complex  extensive 
covalent  mixing  of  metal  3d  and  ligand  2s  and  2p  functions  was 
found.   The  resulting  wave  function  was  used  to  compute  the 
crystal  field  splitting  parameter,  lODq,  and  the  transferred 
hyperfine  interaction   parameters.   Comparison  with  experiment 
led  to  the  conclusion  that  a  Hartree-Fock  description  of  the 
cluster  gave  reasonably  good  charge  densities  and  excitation  energies. 

Unfortunately  these  techniques  are  not  .at  present  applicable 
to  actinlde  complexes.   In  Section  II  below  we  summarize  in  some 
detail  a  recent  variant  of  the  Hartree-Fock  approximation  (the 
self-consistent-f ield-X  multiple  scattering  method)  which  is 
feasable  for  actinlde  complexes  and  in  Section  III  we  report 

some  preliminary  results  for  the  Uranyl  ion(UOp   )  and  Uranium 

7 
Hexafluoride   done  in  collaboration  with  M.  Boring  and  J.  Wood 

of  the  University  of  California  Los  Alamos  Scientific  Laboratory. 


Section  IV  reviews  the  systems  we  propose  for  future  investigation 
and  finally  in  Section  V  we  briefly  speculate  on  the  possibility 
of  using  the  results  obtained  from  the  multiple  scattering  method 
to  guide  the  formulation  of  atomic  pseudopotentials  to  replace 
the  actinide  core  electrons. 


II.   The  Self-Consistent-Field-X  Multiple  Scattering  Method 

o 
The  Hartree-Fock  self-consistent-field  model   is  a  one 

electron  approximation  in  which  the  true  many  electron  wavefunction 

is  approximated  as  an  antisymmetrized  product  of  one  electron  spin 

orbitals,  u..   From  these  spin  orbitals  we  can  set  up  a  charge 

density 

p(l)  =  E^n^u^*(l)u^(l) 

where  the  n.  are  occupation  numbers  with  the  value  zero  or  unity. 
The  expectation  value  of  the  total  energy  then  takes  the  form 

<E>  =  E^n^  /  u^*(l)f^u^(l)dv^  +  \   /  p(l)[  /  p(2)g^2^^2  ^   U(l)]dv^ 

where  f,  is  a  one  electron  operator  consisting  of  the  kinetic 
energy  and  the  potential  energy  of  an  electron  in  the  field  of 
the  nuclei  and  g, „  represents  the  coulomb  repulsion  between  a 
pair  of  electrons.   The  first  term  represents  the  average  of  the 
one  electron  operators.   The  second  term  represents  one  half  the 
Integral  of  the  charge  density  multiplied  by  an  effective  potential 
energy.   This  potential  energy  is  made  up  of  the  coulomb 
interaction,  /  p(2)g-,pdVpj  which  includes  the  self-interaction 
of  the  electron  and  of  the  term  U(l),  a  negative  term  which 
subtracts  the  self-interaction  and  also  takes  account  of  all 
exchange  terms. 

The  one-electron  Schrodinger  equations  used  for  determining 
the  spin  orbitals  u.  are  now  derived  by  applying  the  variation 
principle  to  the  expression  for  the  total  energy. 

4 


[-V^2  +  V^(l)  +  V^^p(l)]u.(l)  =  E.u.(l) 

2 
Here  -V,   is  the  kinetic  energy  and  V  (1)  is  the  electrostatic 

potential  energy,  both  electronic  and  nuclearj  computed  classically. 

The  quantity  V  t,p(l)  is  the  exchange  potential.   Up  to  this  point 

the  formulation  is  exact  within  the  framework  of  the  one  electron 

model. 

The  Hartree-Fock  equations  may  be  solved  exactly  by 
numerical  integration  for  atomic  systems.   For  molecules  they 
may  be  solved  approximately  by  expanding  the  molecular  spin 
orbitals  in  terms  of  a  finite  linear  combination  of  atomic 
orbitals  (Gaussian  or  Slater)  centered  on  the  constituent  atoms. 

At  this  point  two  further  approximations  to  the  Hartree-Fock 
equations  have  been  introduced  which  allow  their  practical 
extension  to  actinide  complexes.   The  first  is  the  statistical 

Q 

approximation  for  the  exchange  potential  V  „„ , 

Xnr 


V^  (1)  =  -6a(3/8TTp(l))^/3 


where  a  is  a  constant  to  be  determined.   The  criterion  used  is 
that  a  should  give  as  good  a  description  as  possible  of  the 
separated  atoms.   Slater  and  Wood   have  shown  that  if  one 
chooses  a  so  that  <E   >  agrees  exactly  with  the  Hartree-Fock 
energy  of  the  same  atom  then  the  X   spin  orbitals  agree  closely 
with  the  Hartree-Fock  spin  orbitals. 

The  second  approximation  introduced  by  Slater,  Johnson,  and 
Connolly    is  to  approximate  the  one  electron  potential  by  the 
so-called  "muffin-tin"  approximation.   That  is,  each  atom  is 

5 


surrounded  by  a  sphere  (Region  I)  and  inside  each  of  these 
spheres  the  exact  potential  is  replaced  by  a  spherical  average. 
In  the  region  between  the  spheres  (Region  II)  the  potential  is 
replaced  by  a  constant.   The  entire  molecule  is  now  surrounded 
by  a  sphere  and  in  this  extra-molecular  region  (Region  III)  the 
exact  potential  is  replaced  by  a  spherical  average. 

In  carrying  through  a  self-consistent-field  calculation  one 
may  start  with  an  assumed  potential  (i.e.,  a  superposition  of 
atomic  potentials),  solve  the  one  electron  equations  for  that 
potential,  finding  spin  orbitals  u.  with  eigenvalues  E..   The 
charge  density  arising  from  the  occupied  spin-orbitals  then 
determines  a  new  potential  used  as  the  starting  point  of  the 
next  iteration. 

The  computational  method  used  for  solving  the  one  electron 

12 
Hartree-Fock  equations  was  formulated  by  Johnson    and  is  known 

as  the  multiple  scattering  method.   Inside  each  atomic  sphere  in 

Region  I  there  is  a  spherically  symmetric  potential.   The  solution 

of  Schrodinger ' s  equation  can  be  expanded  as  a  linear  combination 

of  spherical  harmonics  and  numerically  determined  radial  functions 

depending  on  the  angular  momentum  number  I    and  the  energy  E . . 

The  spherically  symmetric  Region  III  outside  the  molecule  may  be 

treated  in  a  like  manner.   In  the  Region  II  (the  intersphere 

region)  there  is  a  constant  potential.   The  solution  can  be 

taken  to  be  a  sum  of  free-space  scattering  functions  Cp^^oducts 

of  spherical  harmonics  and  the  proper  combinations  of  spherical 

Bessel,  Neumann,  or  Hankel  functions)  which  fall  off  in  the  proper 


exponential  manner  as  one  goes  away  from  an  atom.   The  satisfaction 
of  the  proper  boundary  conditions  on  the  surface  of  the  atomic 
spheres  leads  to  the  determination  of  the  unknown  expansion 

coefficients  as  well  as  the  eigenvalues  E. . 

13 
These  equations  have  been  programmed  by  Johnson  and  Connolly 

and  can  be  solved  very  efficiently.   The  expansion  of  the  wave- 
function  converges  very  rapidly  as  not  many  values  of  angular 
momentum  are  required  for  each  atomic  sphere.   The  final  wave- 
function  has  the  physical  features  of  a  linear  combination  of 
atomic  orbitals.   In  the  traditional  Hartree-Fock  approach  one 
obtains  these  molecular  orbitals  as  a  highly  approximate  solution 
to  the  correct  one-electron  Hamiltonian.   In  the  multiple-scattering 
method  one  solves  a  model  Hamiltonian  exactly  with  great  speed 

and  efficiency.   It  has  recently  been  shown  by  Rosch,  Klemperer, 

14 
and  Johnson    that  approximate  non-muffin-tin  corrections  can  be 

introduced  into  the  method  through  the  use  of  overlapping  atomic 

spheres.   It  should  also  be  pointed  out  that  the  transition 

energies  calculated  in  the  multiple-scattering  method  are  weighted 

averages  of  the  multiplets  arising  from  a  given  atomic  or  molecular 

configuration . 

The  multiple-scattering  method  has  been  applied  with  excellent 

results  to  a  variety  of  molecular  systems.     In  this  laboratory 

we  have  recently  completed  successful  Investigations  of:  (1)  the 

15 
rydberg  states  of  water   ;  (2)  the  electronic  spectra  of  the 

16 
solvated  electron   ;  (3)  the  electronic  spectra  of  various 

17 
halogen  substituted  permanganates  and  chromates   ;  and  (4)  a 

2-  18 
study  of  metal-metal  bonding  in  Re2Clo 


III.   Preliminary  Results 

l4 
Over  a  decade  ago  McGlynn  and  Smith    interpreted  the 

electronic  structure  of  the  actinyl  ions  based  upon  a  molecular 

orbital  picture.   They  postulated  the  ground  state  of  the  linear 

^         ion    to   be   of      Z 
2  g 


UOp         ion    to   be   of      Z      symmetry   with    configuration. 


,      2  ,      2  ,      i|  ,      4 

la„    la ^   iTT,,    liT^ 

u        g        u        g 

The  molecular  orbitals  were  constructed  from  2p  orbitals  on  each 

oxygen  and  5fj  6d,  7s,  and  7p  orbitals  centered  on  the  uranium 

atom.   In  the  Russell-Saunders  (Hund's  case  a)  coupling  scheme 

the  first  excited  state  arises  from  either  the  configuration 

■^l      ->  ITT  2  (|)  -'-  (-^'^A  ,  -^'^r  )  which  is  orbitally  forbidden  or  the 
g     g  ^u      u'      u  ^ 

configuration   Z   ->  lir  --^  6  (  '  IT  ,   '  $  )  the  Z  -*■  IT  component  of 
g     g   u     u      u 

which  is  allowed.   Upon  analysis  of  the  U0„    absorption  and 
emission  spectra  in  solution  McGlynn  and  Smith  assign  the  first 

three  bands  to  the  components  of  a   n   state. 

20 
Newman    has  recently  interpreted  the  near  visible  absorption 

spectra  of  UOp    in  the  cesium  uranyl  nitrate  crystal.   He 

emphasizes  the  importance  of  spin-orbit  coupling  and  suggests  that 

jj -coupling  (Hund's  case  c)  is  dominant  in  the  uranyl  ion.   In 

addition  Newman  postulates  a  different  orbital  ordering  in  the 

ground  state  based  on  U-0  overlap  integrals  computed  with 

4   2    4    2  ^ 
numerical  Hartree-Fock  function  (Itt   la   Itt   la   ). 

21 
Most  recently,  Gorller-Walrand  and  de  Jaegere    reported  a 

systematic  study  of  a  large  number  of  uranyl  complexes.   The 

visible  absorption  spectra  of  these  complexes  could  be  classified 

8 


into  six  distinct  categories  differing  in  the  number  and  relative 
intensity  of  the  electronic  transitions  involved.   To  each  of 
these  spectral  categories  could  be  associated  a  given  geometry. 

By  analyzing  these  results  by  means  of  group  theory  and  simple 

22 
molecular  orbital  theory,  Gorller-Walrand  and  Vanquickenborne 

have  shown  that  the  band  system  situated  at  20000-22000  cm~"^  is 

only  consistent  with  a   S   -*-   $   transition  (a   (j)   )  ,  while  the 

band  system  at  22500-27000  cm""""  is  consistent  with  a  "'"I   -^  "^A 

transition  (a  ^^    ^  •      ^"^   addition  they  have  shown  that  it  is 

possible  to  use  the  small  equatorial  perturbation  due  to  the 

ligands  as  a  physical  probe  to  examine  the  coupling  situation 

23 
in  the  excited  uranyl  $   and  A   states.     They  conclude  that 

only  a  Russell-Saunders  coupling  scheme  allows  a  satisfactory 

explanation  of  the  observed  splitting  and  intensity. 

7 
In  our  computations  on  the  uranyl  ion   a  starting  potential 

for  the  multiple  scattering  method  was  generated  by  superimposing 

atomic  charge  densities.   For  the  U   ion  a  relativistic  atomic 

24 
(Dirac-Slater)  charge  density    was  superimposed  with  a  non- 
10 
relativistic  oxygen  atomic  (V   )  charge  density.     The  valence 

electrons  for  the  U    ion  were  considered  to  be  6s  6p  5f  6d  7s 

while  the  remaining  electrons  of  the  core  were  held  fixed.   In 

p 
the  case  of  the  oxygen  atoms  only  the  Is   electrons  were  considered 

part  of  the  core.   Once  a  starting  potential  was  obtained  by 

superimposed  atomic  potentials  the  calculations  were  taken  to 

self-consistency.   We  summarize  below  the  results  for  linear  UO^ 

at  R  =  1.73  A  (3.269  a.u. ) . 


First,  the  four  highest  occupied  molecular  orbitals 

2j  2   4   2 
fell  in  the  order  ii   o   tt   a   .   These  orbitals  show  heavy 

mixing  between  the  uranium  5f  and  6d  orbitals  and  the  oxygen 

2p  orbitals  with  little  participation  of  the  uranium  7s  or  7p 

orbitals  or  the  oxygen  2s  orbital.   The  covalency  is  thus 

surprisingly  large.   The  total  number  of  electrons  in  each  region 

for  these  four  molecular  orbitals  is  as  follows:   uranium  sphere 

(total  3.66,  .02s,  .Ip,  .58d,  2.96f);  in  each  oxygen  sphere 

(total  2.06,  2.06p);  outer  sphere  (total  .36,  .04s,  .04p,  .l8d,  .lOf) 

intersphere  region  (total  3.86).   Although  it  is  difficult  to 

assign  the  outersphere  and  intersphere  charge  to  a  particular  atom 

most  of  the  charge  probably  belongs  to  the  oxygens.   This  is  due 

to  the  fact  that  the  atomic  spheres  were  not  allowed  to  overlap 

resulting  in  a  far  larger  uranium  (2.33^4  a.u.)  than  oxygen 

(0.93^6  a.u.)  atomic  sphere  radius.   The  uranium  radius  was  set 

by  requiring  it  to  contain  the  entire  core  charge  density.   The 

calculations  thus  show,  depending  upon  how  the  charge  is 

partitioned,  a  slight  ligand  to  metal  charge  transfer. 

We  have  also  examined  the  electronic  transition  energies  for 

the  excitation  -^T.      ^   -^'^^    ( a  "^  (})  "^  )  which  falls  at  l6,457  cm""'-  and 
g       g  u   u 

the  transition  ^Z   ^  "^'^n  ,  ^'^f  (^  ^5  ^)  which  falls  at  219^2  cm"^ . 
g       u      u  g  u 

It  must  be  remembered  that  in  the  multiple-scattering  method  the 
transition  energies  represent  weighted  averages  of  the  multiplets 
arising  from  a  given  molecular  configuration.   The  results  are 

thus  suggestive  but  certainly  not  conclusive  that  the  interpretation 

21-23 
of  Gorller-Walrand  et.  al.       are  correct. 


10 


We  report  briefly  on  some  preliminary  results  on  the  charge 
distribution  In  UPg  and  on  the  optical  spectra  of  UFg~  •   A 

detailed  review  of  the  absorption  spectra  of  actlnlde  compounds 

PR  1 

has  been  published  recently  by  Ryan.     In  addition  5f  transi- 
tions. In  particular  pentavalent  uranium  have  been  reviewed  by 

Selbln  et.  al . ,   and  Rlesfeld  and  Crosby.    The  5f   configuration 

2 
gives  rise  in  octahedral  symmetry  to  a  Ap   ground  state  and  a 

2        2 
Tp   and   T-,   excited  states.   In  the  intermediate  coupling  scheme 

2         2 
the  ground  state  is   r„  while   T„   is  split  by  spin  orbit  coupling 

into  ^Fg  and  ^I'   and  ^T-^^  is  split  into   T'g  and  ^Tg. 

Experimentally  we  would  thus  predict  at  least  four  low  energy 

transitions  to  be  observed,  although  Jahn-Teller  splitting  of 

29 
the  fn  and  T'o  states  is  possible.   Rlesfeld  and  Crosby 

assign  the  spectrum  as  follows:   Ap  ^  Tp  (rg,  f  )  at  ^587  cm" 

and  6928  cm"^;  A„   ->  T,  (T'p,  T  r)    at  12705  cm~^  and  1^1245  cm"""- . 

2U      lU     CS     D 

The  calculation  we  performed  on  uranium  hexafluorlde  was 
carried  out  in  an  essentially  equivalent  manner  to  that  reported 
above  for  the  uranyl  ion.   All  electrons  below  the  6s  and  6p 
level  of  uranium  were  treated  as  part  of  the  core  along  with 
the  Is  electrons  of  the  fluorine  atoms.   The  covalent  mixing  of 
the  metal  5f  orbltals  with  the  llgand  2p  orbltals  was  again 
surprisingly  large.   The  distribution  of  the  56  valence  electrons 
was  as  follows:   uranium  sphere  (total  10.72,  1.96s,  5.19p,  .7d, 
2.87f);  In  each  fluorine  sphere  (total  5-55,  1.62s,  3.93p); 
outersphere  (total  .75,  .05s,  .15p,  .27d,  .28f);  intersphere 
region  (total  11.25).   The  2 . 87f  electrons  found  in  the  uranium 
sphere  may  be  contrasted  to  the  3-7f  electrons  found  by  Hill  and 


11 


4 
Kmetko   as  a  result  of  their  band  structure  calculation  for 

metallic  uranium.   We  have  also  performed  a  series  of 

computations  upon  UF^r"  and  computed  the  one  electron  excitation 

energies  A^^  ^  T^^C^SSS  cm"^)  and  A^^  ->  T^^  (17005  cm"^). 

These  are  seen  to  be  in  fair  agreement  with  the  experimental 

results.   It  should  be  remembered  that  within  the  multiple 

scattering  formulation  one  electron  transitions  of  this  type  are 

computed  exactly;  i.e.,  there  is  no   problem  arising  from  the 

separation  of  multiplets. 
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IV.   Proposed  Investigations 

The  results  of  the  calculations  reported  above  strongly 
suggest  the  possibility  of  obtaining  both  qualitative  and 
quantitative  information  about  the  electronic  structure  and 
bonding  of  actinide  complexes  by  means  of  the  self-consistent- 
field  multiple  scattering  model.   The  recent  work    on  over- 
lapping atomic  spheres  gives  hope  of  partially  sidestepping  the 
effects  of  the  muffin-tin  approximation  and  thus  reducing  the 
number  of  electrons  assigned  to  the  intersphere  region.   One 
should  also  investigate  the  changes  that  might  be  expected  if 
relativistic  corrections  were  included  for  the  valence  electrons 
as  well  as  for  the  frozen  core.   It  should  be  possible  to  put 
the  relativistic  corrections  directly  into  the  Hamiltonian  as  is 
done  in  relativistic  band  structure  calculations.     We  propose 
here  as  a  first  step  in  that  direction  to  compute  the  spin-orbit 
splitting  of  the  upper  few  molecular  orbitals  as  a  perturbation 
correction  in  a  manner  similar  to  that  used  by  Wood    in  atomic 
X   computations. 

We  now  summarize  briefly  some  of  the  systems  that  seem 
fruitful  for  future  investigation. 

A.   Actinyl  ions 

We  propose  to  complete  the  series  of  calculations  begun 
with  U0„    to  include  NpO„   ,  PuOp   ,  and  AmOp    whose  spectra 
are  reasonably  well  characterized.   These  systems  are  small 
enough  to  allow  experimentation  with  the  method  in  order  to 
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develop  inproved  computational  techniques.   This  is  particularly 
true  for  NpO„    where  the  configuration  5f   allovrs  simple 
interpretation  of  the  spectra  in  the  multiple-scattering  model. 

B.   Actinide  Hexahalides  and  Oxides 

The  reviews  of  Ryan    (halides),  Selbin  and  Ortega    (U   ), 
and  the  recent  work  of  Hecht,  Lewis,  and  Eastman   (oxides)  indicate 
the  extraordinary  range  of  experimental  information  available  on 
the  electronic  structure  of  actinide  complexes  and  in  particular 
those  with  the  configuration  5f  •   We  propose  to  continue  our 
-investigation  of  UFg  and  UF^~  and  extend  the  computation  to 
other  halide  ligands .   In  addition  we  plan  to  investigate  NpF^- , 
the  5f^  oxides  of  uranium  (UO^,  UO^"^,  UO  ~ ,  UO^^^",  and  \]^Q^^~), 
and  several  compounds  of  plutonium.   It  is  proposed  to  perform 
spin-polarized  computations  in  order  to  compute  the  ligand 
transferred  hyperfine  parameters. 


C.   Metal  Clusters 

The  actinide  metals  exhibit  a  wide  variety  of  interesting 
magnetic,  optical,  and  superconducting  properties.   In  recent 
years  a  number  of  band  structure  calculations  have  been  performed 

in  order  to  gain  better  theoretical  understanding  into  their 

Q  4  [5 
origin  '  '  .   Unfortunately  the  metals  have  complex  crystal 

32 
structures  with  many  atoms  per  unit  cell  and  low  symmetry 

Thus  recent  computations  on  uranium  have  considered  only  such 

phases  as  y-uranium  with  body  centered  cubic  (bcc)  structure. 
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We  propose  to  perform  a  self-consistent-field  multiple-scattering 

calculation  on  a  nine  atom  bcc  uranium  cluster  and  compare  the 

4 
results  to  the  band  structure  calculations  of  Kmetko  and  Hill 

in  order  to  calibrate  our  results.   It  will  then  be  possible  to 

perform  computations  on  the  five  atom  cluster  of  a-uranium 

(Cp   symmetry).   The  extension  of  this  technique  to  the  atomic 

clusters  arising  from  the  different  phases  of  neptunium  and 

Plutonium  would  be  straightforward  and  should  yield  interesting 

information  about  the  participation  of  5f  electrons  in  different 

bonding  situations. 
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V.   Atomic  Pseudopotentials 

The  multiple-scattered-wave  formalism  developed  by  Slater 
and  Johnson  is,  as  discussed  above,  the  only  practicable  method 
which  may  be  applied  to  the  study  of  the  electronic  structure  of 
actinide  complexes.   Its  great  strength  is  the  ease  and  economy 
with  which  it  may  be  applied  to  systems  containing  heavy  metals. 
The  predictions  of  charge  distributions  and  photo-electron 
spectra  have  been  encouraging.   The  prediction  of  optical 
properties  is  more  difficult  as  the  method  is  only  capable  of 
separating  multiplets  in  situations  of  low  symmetry.   The  method 
seems  least  appropriate  when  discussing  stability  or  conformation 
of  molecular  complexes.   The  necessity  to  compute  total  energy 
arises  in  the  prediction  of  vibrational  force  constants  and 
anharmonicity  constants  important  in  the  discussion  of  isotope 
separation,  for  example,  that  of  UFg. 

An  alternate  technique  to  the  multiple-scattering  method 
would  be  to  replace  the  core  electrons  of  the  metal  with  an 
atomic  pseudopotential   .   We  are  in  this  laboratory  at  present 
engaged  in  implementing  a  set  of  atomic  pseudopotential  programs 
based  on  the  work  of  Goddard,  Kahn,  and  Melius""  .   If  we  can 
successfully  apply  this  meth^od  to  transition  metal  complexes  it 
should  be  possible  to  extend  the  formulation  to  the  actinides. 
However,  unlike  the  transition  metal  complexes  which  are 
reasonable  well  understood,  it  will  be  crucial  to  have  the  results 
of  th^e  multiple-scattering  computations  to  guide  the  application 
of  pseudopotential  methods.   We  propose  to  follow  this  promising 
line  or  research. 
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